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Temperature crossover of decoherence rates 
in chaotic and regular bath dynamics 
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The effect of chaotic bath dynamics on the decoherence of a quantum system is examined for the 
vibrational degrees of freedom of a diatomic molecule in a realistic, constant temperature collisional 
bath. As an example, the specific case of U in liquid xenon is examined as a function of temperature, 
and the results compared with an integrable xenon bath. A crossover in behavior is found: the inte- 
grable bath induces more decoherence at low bath temperatures than does the chaotic bath, whereas 
the opposite is the case at the higher bath temperatures. These results, verifying a conjecture due 
to Wilkie, shed light on the differing views of the effect of chaotic dynamics on system decoherence. 
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I. INTRODUCTION 



Decoherence, and the control of decoherence, is a cen- 
tral problem in modern quantum physics. In particu- 
lar, "quantum technologies" .such as quantum cryptog- 
raphy, quantum computing [l[ and quantum control [2j 
rely upon the maintenance of quantum effects over sig- 
nificant periods of time. As such, decoherence serves as a 
primary obstacle to progress in the experimental imple- 
mentation of a number of these quantum based scenarios. 
Of particular interest is the nature and rates of decoher- 
ence in systems in the condensed phase. 

In this regard, it has been long argued |3l4l4]| that 
traditional uncoupled oscillator or standard spin bath 
models (e.g., the spin-boson [15|, [lf|, boson-boson J171 - 
Il9| . spin-spin |20j models) are inadequate to describe 
dynamics in condensed phases since they lack intra- 
cnvironmental coupling. Such intra-environmental cou- 
pling can display different types of behavior, including 
chaotic dynamics, and therefore its effect on decoherence 
can be significantly different from traditional models. 

For example, in such uncoupled bath cases, when the 
system perturbs the bath, it cannot relax internally; en- 
ergy must flow through the subsystem in order for the 
bath to return to equilibrium. This causes the system 
energy to increase initially even when relaxation is ex- 
pected [21(. Second, the system becomes strongly en- 
tangled with the bath as a result of this energy flow and 
hence the system decoheres more strongly than it should. 
Finally, the equilibrium state that the bath reaches may 
not be the expected canonical state [221 ] . 

Applications of uncoupled oscillators to model con- 
densed phase environments may be especially problem- 
atic. For example, anharmonic corrections are known to 
be important in the study of phonons in Si [23| and are 
essential for the explanation of heat transport. Both the 
vibrational dynamics of Si [2J] and its electronic struc- 
ture [25| are believed to be chaotic. Dynamics of a col- 



loidal particle in water have also been shown to be chaotic 

The effect of the structured environment on some is- 
sues in solid state has been examined. For example, it has 
recently been shown, for the central spin model [27H33J . 
that the role of structured environments on solid-state 
(ferromagnetic phase) implementations is important. In 
particular, the dynamical regime of the bath has been ob- 
served [28| to determine the efficiency of the decoherence 
process. For example, in a perturbative regime, deco- 
herence is stronger in the integrable limit; on the other 
hand, in the strong coupling regime the chaotic limit is 
more efficient. Also, the two-spin system decoherence has 
been found to exhibit different behavior depending on the 
characteristics of the coupling with the environment, as 
well as on the internal dynamics and initial state of the 
environment. 

Clearly, studies on more realistic models than the un- 
coupled oscillator or spin models are needed for con- 
densed phase environments. Efforts to generalize the os- 
cillator bath model in such cases are very preliminary. 
When intra-bath coupling is added to the boson or spin 
bath, general analytic solutions are unavailable and ex- 
act solutions can only be obtained computationally for 
very small baths. In the case of oscillator bath models, 
the bath cannot consist of more than three or four os- 
cillators. For spin baths results have been reported for 
approximately 20 spins. 

The most common question addressed in such studies 
is whether intra-bath coupling increases or decreases the 
decoherence of the embedded subsystem. Surprisingly, 
this apparently simple question has generated consider- 
able controversy. It was first conjectured by Zurck that 
decoherence should be greater for chaotic baths [3| . This 
was quantitatively verified in a study wherein a single 
harmonic oscillator subsystem interacts with a bath con- 
sisting of a single chaotic oscillator [5]. Unfortunately, 
the relevance of this result to the usual paradigm of a 
small system interacting with a large environment is un- 



clear. Furthermore, Alicki has argued [7(], by contrast, 
that the decoherence rate in the limit of pure decoher- 
ence (i.e. in the absence of dissipation) will be greater 
for an integrable bath than for a chaotic bath. This is 
because the energetically available bath states in the in- 
tegrable case can be highly degenerate, whereas only one 
state is available in the chaotic bath at a given energy 
due to level repulsion. Thus if a chaotic system-bath 
state is microcanonical, the wavefunction of the system 
plus bath will be a simple product state, since there is 
only one energetically available bath state with which the 
system can couple. Accordingly, there should be greater 
decoherence in the integrable case. This viewpoint can 
also be supported with semiclassical arguments [a] . It can 
be shown that the square of the off-diagonal matrix ele- 
ments of the system-bath coupling operator scale as h 1 ^" 1 
for a chaotic bath [34| . The off-diagonal coupling matrix 
elements thus vanish in the thermodynamic limit. At 
low temperatures the diagonal matrix elements change 
slowly with energy [6( and so the subsystem dynamics is 
shifted but not strongly decohered. By contrast, selection 
rules for integrable systems guarantee large off-diagonal 
matrix elements which cause strong decoherence. These 
conclusions were verified numerically for a low tempera- 
ture spin-bath @. 

Thus, there are two well-defined and apparently con- 
tradictory positions on the issue of whether a chaotic 
bath may increase or decrease decoherence. Numerous 
low temperature spin-bath studies |4.la.l8l- ll2l . [l4l | support 
Alicki's [7| predictions that decoherence should be greater 
for integrable baths. However, there are spin-bath stud- 
ies that draw the opposite conclusion [13J and hence sup- 
port Zurek's conjecture that chaotic baths cause greater 
decoherence. 

Based on these results Wilkie has speculated [35| that 
some sort of transition occurs with increasing temper- 
ature; a chaotic bath could cause less decoherence at 
low temperature and greater decoherence at high tem- 
perature. This conjecture would be consistent with the 
spin-bath results 0, |g, I8l4l4j and would not be in direct 
conflict with the oscillator calculation [5]. However, the 
existence of such a transition is difficult to verify in exact 
spin-bath or oscillator-bath calculations. High tempera- 
ture calculations for a bath of ten spins [6y did not show 
such a transition and calculations for larger spin-baths 
could not be carried out at high temperature. 

An alternative approach would be to explore the possi- 
bility of such a transition using an approximation scheme. 
This approach is the focus of this paper. Recently it has 
been shown that quantum decoherence can be accurately 
computed using classical dynamics simulations based on 
the quantum Wigner function [3a . [37j . In this Wigner 
approach, regions of phase space where the Wigner func- 
tion of the initial state takes negative values are Monte 
Carlo sampled using the absolute value and the resulting 
classical trajectories carry a negative sign as a weighting 
factor. The resulting approximation can be very accu- 
rate, and in this paper we employ this approach to ex- 



amine the decoherence of a superposition of vibrational 
states of I2 in liquid Xe baths comprising 512 atoms. 

Our key observations are: (a) we observe less deco- 
herence of vibrational superposition states of I2 at low 
temperatures for liquid Xe than for its ideal gas coun- 
terpart obtained via simulations without Xe-Xe interac- 
tions, but (b) as the temperature is increased, a tran- 
sition is observed and liquid Xe becomes the stronger 
source of decoherence. Thus, we show the existence of 
the two regimes in a physically realistic model. Note, as 
an immediate application, that the decoherence of a vi- 
brational superposition state is a significant impediment 
to coherent control via pump-dump scenarios [2j. Hence, 
understanding conditions responsible for decoherence in 
such systems is important. Indeed, this was the original 
motivation for examining this particular system. 

This paper is organized as follows. The model consid- 
ered, as well as the details of the numerical simulations, 
are discussed in Sec. ITT] Section Hill reports the numerical 
results of the simulations for two different initial states at 
three temperatures and qualitative explanations for the 
observed behavior are proposed. Finally, Sec. lIVI contains 
a summary of conclusions. 



II. MODEL 

A. Hamiltonian 

Consider the decoherence rates for different superpo- 
sitions of vibrational states of I2 coupled to a bath of 
Xe atoms. The subsystem of interest is the vibrational 
degree of freedom of the diatomic, and the environment 
comprises the translational degrees of freedom of the I2 
and of the Xe atoms. The Hamiltonian describing the 
full system [38j can be written as a standard system- 
plus-environmcnt Hamiltonian, as follows: 



where 
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Here, Eqs. (j2"aj) and (l2b|) describe the independent evo- 
lution of the I2 vibrational degree of freedom (q) and the 
bath dynamics, respectively, while Eq. (j2"c) accounts for 
their interaction. The isolated I2 is described by a Morse 
oscillator, 



V = D 1 



3 -/3(9-<3o) 



(3) 



with D = 1.2547 x 10 4 cm- 1 , (3 = 1.8576 A" 1 , q = 0, 
and p\ 2 = mi 2 /4, with ra\ 2 being the I 2 molecule mass 
and p\ 2 being its reduced mass. The degrees of freedom 
of the environment include the translational degree of 
freedom of the 1 2 (its center of mass, with mass mi 2 ), as 
well as the collection of N Xe atoms. The interaction 
between Xe pairs is described by the Xe-Xe interaction 
potential <pXe-Xe{fij) j where r^ is the distance between 
the ith and jth Xe atoms. This interaction is modeled 
by a realistic pairwise Lennard- Jones potential: 



OXo-Xc 



(nj) = 4e Xc -j 



0"Xc-Xc 



12 



CXo-Xc 



( 4 ) 
where exc-Xc = : 154.00 cm 1 is the well-depth of the 

potential and oxc-Xc = 3.930 A is related to the posi- 
tion of the minimum of the well [V(r m i n ) = — exc-Xc], 
r m j n = 2 1 ' 6 crxc-Xc- Here, exo-Xe gives an estimate of 
the intensity of the interaction between two Xe atoms 
and (Txc-Xc is the effective diameter of the Xe atoms. 
Note that under these conditions and at the densities 
and temperatures considered in this work Xe is a liquid. 

The coupling between the I2 vibration and the bath is 
described by the interaction Hamiltonian [see Eq. (|2Tc|) ]. 
If vibration and translation were not coupled in the I2 , 
this term would just account for the interaction between 
the I2 with the Xe atoms and, therefore, would look 
like Eq. ©, with e Xo -Xo and oxc-Xc replaced by ei 2 -x c 
and <7i 2 _xc: and r^ replacing r oi between the I2 and 
the iih Xe atom. In such a case, £i 2 _xc an d ci 2 _x c 
can be taken as the average of the corresponding Xe- 
Xe and I2-I2 interactions, i.e., e = ^/exo-Xcei 2 -i2 an d 



cr, 



(0) 
I 2 -Xe 



(crxe-Xe+ui 2 _i 2 )/2 (here, ei 2 -i 2 = 382.27cm 



and cti 2 -i 2 = 4.982 A denote the well-depth and posi- 
tion of the minimum for the corresponding I2-I2 pair- 
wise Lennard- Jones potential function). Here, however, 
since the diatomic is "breathing" while it vibrates, a 
plays the role of an effective radius given by cti 2 -xc = 



(oxe 



-Xo 



0I2-I2 



aq)/2 



r (0) 



'i 2 -Xc 



aq/2 [38|. Since 
diatoms expand and contract in one direction, a < 1. 
This model, termed an effective breathing sphere, models 
the interaction of I2 with the surrounding environment 
through the interaction potential: 



(f>i 2 -x c (r Qll q) = 4e 



I2-X0 



Cl 2 -Xo 

r i 



12 / x 6 

CI2-X0 



roi 



(5) 

a potential that is known [38j to be quantitatively reliable 
for this system. 

To obtain the analog of an uncoupled oscillator bath 
for comparison purposes, we simply ignore the Xe-Xe in- 
teractions, setting the third term in Eq. (I2bl) to zero. 
Collisional interactions between the Xe atoms of the bath 
are thus removed, although collisions with the I 2 are re- 
tained, a model referred to below as the "Xe ideal gas" . 



B. Dynamics 

We consider initial states consisting of a thermally 
equilibrated Xe bath within which is embedded an io- 
dine molecule in a superposition of vibrational states. 
Conceptually, such a superposition could be prepared by 
laser excitation from the ground vibrational state, where 
a multiphoton path can be utilized to overcome any se- 
lection rule issues. The subsequent dynamics calcula- 
tions are done by sampling the Wigner distributions cor- 
responding to the initial superpositions of Morse eigen- 
states 39] with classical trajectories [4y|, followed by 
standard Molecular Dynamics (MD) techniques [4J| us- 
ing the velocity- Verlet algorithm [42[ to propagate the 
trajectories. The result is the time-evolving density 
p(q,p,r ,p ,{ri,p t }f =1 ,t) for the full system + bath. A 
total number of 2 x 10 6 trajectories was considered. In 
all reported simulations the total number of particles (I2 
plus Xe atoms) is 512, which was found to be adequate 
to converge the calculation of the purity, used as a mea- 
sure of the system coherence. Finally, a density (I2 in 
Xe) p* — 3053 g/cm 3 was used in all the calculations 
to fix the size of the MD cell. Conversion factors and 
parameters used in the simulation are provided in the 
Appendix. 



III. COMPUTATIONAL RESULTS AND 
DISCUSSION 

As a measure of decoherence, we compute the purity 
X of the I2 dynamics [43[ , defined as 



X = Tr[Ps(*)] = / p 2 s {q,p,t)dqdp, 



(6) 



where 



Ps(q,P,t)= / p{q,P,rQ,p ,{r i ,p i \ i=1 ,t)dr dp \i i=1 dr. l dp, 

(7) 
is the reduced density associated with the subsystem of 
interest, which here is the I2 vibrational degree of free- 
dom. Since the initial state (a vibrational superposition) 
is described by a wavefunction, the purity is initially 
unity, but decays with time as a consequence of the en- 
tanglement of the system and bath degrees of freedom. 
A small amount of decay in \ is a l s ° observed as a func- 
tion of time for the isolated diatomic propagated in the 
absence of the bath. This decay [see, e.g., Fig. [Tj is a 
measure of the computational accuracy and is found to 
be very small over the relevant 5 ps time scale. 

Figure HJa) shows the purity as a function of time for 
the two cases of Xe liquid and Xe ideal gas where the 
initial vibrational degree of freedom of the I2 is in an 
equal superposition of the ground and second excited vi- 
brational state. At 177.36 K both liquid Xe and ideal gas 
Xe are seen to cause decoherence of I2 on a picosecond 
time-scale. In this case it is apparent that the liquid Xe 
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FIG. 1: Dependence on temperature of the purity as a func- 
tion of time for the I2 initially in a superposition of the ground 
and second-excited states: (a) T = 177.36 K, (b) T = 221.7 K 
and (c) T = 554.25 K. In all graphs, the black solid line in- 
dicates the Xe-Xe coupling is active, the red dashed line cor- 
responds to the "ideal gas" bath (no intra-bath interactions), 
and the blue dotted line corresponds to the isolated I2. 
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FIG. 2: Dependence on temperature of the Lyapunov expo- 
nent as a function of time for the I2 initially in a superposition 
of the ground and second-excited states: (a) T — 177.36 K, 
(b) T = 221.7 K and (c) T = 554.25 K. In all graphs, the black 
solid line indicates the liquid Xe case and the red dashed line 
indicates the ideal gas Xe results. Note how A(t) is similar in 
all three cases for the ideal gas Xe bath case. 



bath causes substantially less decoherence than does the 
ideal gas Xe bath. 

To verify that the liquid Xe dynamics is indeed chaotic 
we calculated the Lyapunov exponent, denoted by A(i), 
as the distance between two nearby Xe initial condi- 
tions [44j| . which is known to grow exponentially for a 
chaotic distribution dynamics and sub-exponentially for 
integrable dynamics. Examination of the Lyapunov ex- 
ponent A(i) as a function of time for the two cases [see 



Fig. [H[a)] clearly shows exponential growth for liquid Xe 
and sub-exponential growth for ideal gas Xe. Hence, the 
liquid Xe is chaotic, whereas the "ideal gas" case is not. 
The results shown in Figure Q] support the Alicki conjec- 
ture that bath chaos tends to suppress decoherence at 
low temperatures. 

Similar results to those above are obtained at a temper- 
ature of 221.7 K. Figures [T{b) and[5Jb) show the purity 
and Lyapunov exponents, respectively, plotted against 
time for the two cases for the same 0-2 superposition. 



The liquid Xe shows the expected exponentially increas- 
ing dependence of A(i) with time, showing that the bath 
is chaotic. Indeed, the slope of the logarithm of A(i) has 
increased, suggesting that the liquid is even more chaotic 
than at the lower temperature. Again, the liquid Xe case 
is found to induce less decoherence than does the ideal 
gas Xe, but the difference between them is now less pro- 
nounced. 

The situation at high temperature, however, is quite 
different. As seen in Fig. [TJc), at 554.25 K liquid Xe 
bath causes greater decoherence than does the ideal gas 
Xe. The Lyapunov exponent displayed in Fig. HJc) con- 
firms that liquid Xe bath is chaotic, with an even larger 
Lyapunov exponent [as manifest in the slope of A(t) ver- 
sus t]. Thus, a transition has occurred from a low tem- 
perature regime where the bath chaotic bath induces less 
decoherence than does the non-collisional bath, to a high 
temperature regime where the reverse is the case. 

It is evident from Fig.|3]that the decoherence dynamics 
for the integrable "ideal gas" case is relatively unchanged 
as a function of temperature. Hence, it is the chaotic 
case that goes from weak to strong decoherence as the 
temperature increases, and is responsible for the cross- 
over behavior in the decoherence of regular versus chaotic 
baths. 

To verify that this transition occurs for other initial 
conditions we also show simulations for a second coher- 
ent vibrational superposition state. Figures[3ja)-(c) show 
the purity as a function of time for three different tem- 
peratures (177.36, 221.7, and 554.25 K) calculated for a 
coherent superposition of the fifth and eighth vibrational 
states of the I2. The decay of the purity is faster for this 
second initial state than it is for the first. Here the dy- 
namics appears to occur on two time-scales, a fast initial 
decay followed by a slower falloff. The fast initial decay 
is common to both the liquid and ideal gas simulations. 
Otherwise the plots are qualitatively similar to those for 
the first initial state examined above. Significantly, the 
cross-over at high temperatures is again apparent: Deco- 
herence of the liquid is smaller than that of the ideal gas 
at 177.36 and 221.7 K, but larger at the highest temper- 
ature 554.25 K. 



IV. CONCLUSIONS 

We have explicitly compared decoherence in the same 
system immersed in a liquid Xe bath, a well known 
paradigmatic system. In doing so we have confirmed that 
subsystems interacting with condensed phase environ- 
ments cannot be analyzed using approaches that neglect 
intra-environmental coupling within the bath. Further, 
we demonstrated that this coupling leads to a cross-over 
in the decoherence dynamics as a function of tempera- 
ture. 

The transition described in the previous section has 
not been previously observed computationally. It is clear 
that the integrable bath case, in the liquid domain, is 




1.0- 
0.8- 






^ 0.6- 




0.4- 


j\ 


0.2- 


1 ^ 


0.0- 


1 I 1 I 1 I 1 I 1 I 1 



12 3 4 5 6 
t(ps) 

FIG. 3: Dependence on temperature of the purity as a func- 
tion of time for the I2 initially in a superposition of the fifth- 
and eighth-excited states: (a) T = 177.36 K, (b) T = 221.7 K 
and (c) T = 554.25 K. In all graphs, the black solid line 
indicates the Xe-Xe coupling is active, the red dashed line 
indicates the off coupling situation, and the blue dotted line 
indicates the isolated I2. 



relatively insensitive to temperature whereas the chaotic 
shows increasing decoherence with increasing tempera- 
ture. Since the chaotic bath case displays weaker coher- 
ence than the integrable bath case at low temperature, 
the increasing decoherence of the chaotic case with tem- 
perature results in a cross-over of behavior. The low 
decoherence of the chaotic case at low temperatures is in 
accord with the analysis in terms of spectral properties of 
the bath. What is remarkable to note is that the Wigner 



method is capable of properly displaying this behavior. 

The results suggest a more detailed analysis of the ori- 
gins of the difference between the chaotic and integrable 
baths in different density regions and the interrelation- 
ship between the rates of decoherence and the Lyapunov 
exponent of the chaotic bath would be of interest. Work 
on such systems is planned. 
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3.32 ps) and where the magnitudes with asterisk denote 
the re-scaled magnitudes. Thus, for example, a density 
p* = 0.85 will corresponds to p = 3.053 g/cm 3 and a 
temperature T* = 1.26 to T = 280 K, while Planck 's 

constant will become h* = fi/ Jinxc-Xcffxc-Xc^c-Xc — 
0.010388. This scaling leads to a system of dimensionless 
equations of motion, which are solved by means of the 
standard velocity- Verlet method in the case of both the 
system and the environment. Within this scheme, quan- 
tum features are taken into account initially in terms 
of the classical Wigner method, i.e., by carrying out a 
Monte Carlo sampling based on the Wigner distribution 
of the initial state of the I2 vibrational state. 



Computational aspects 

The MD simulations were carried out using well-known 
standard procedures [411 ]. All the parameters involved 
were re-scaled taking into account the parameters asso- 
ciated with the solvent particles (here, the Xe atoms). 
That is: 



interparticle distance: 


r* = r/axc-Xc, 


time: 


t* = i]t, 


frequency: 


uj* = u/r), 


density: 


P* = ^Xc-XcP' 


temperature: 


T* = k B T/e Xc -x c 

1 In 



with rj = Jexc-Xc/mx c cr' Xc _ Xc 

(i.e., 1 MD time unit is equivalent approximately to 



The constants of solute and solvent are those pre- 
viously obtained for the Xc + I2 system [38|, i.e., 
ei 2 _i 2 //c B = 550 K, cti 2 _i 2 = 4.982 A, exc-Xc/fcs = 
221.7 K and er Xc -Xc = 3.930 A (k B = 1.3806505 x 
10~ 23 J/K being Boltzmann's constant). The solute and 
solvent masses are, respectively, too = 4.22 x 10~ 22 g and 
to s = 2.18 x 10~ 22 g. Now, taking into account the re- 



Xc-Xc 



1 and e 



Xc-Xc 



1 



1.936, a\ = 1.268 and 



scaling, we will find rn Xc = 1, a 
for the Xe atoms, while % 
6j = 2.481. Regarding the I2 vibrational degree of free- 
dom, we will find that D* = 81.4208, j3* = 7.30037 and 
/ij = 0.4839. The vibrational frequency corresponding 

to the I 2 is wi 2 = y/2(3 2 D/p h = 4.0451 x 10 13 s" 1 , which 
in MD reduced units becomes w T * = 134.16 in reduced 
MD units (the period, in these MD units, is therefore 
t = 2tt/w ~ 0.0468). 
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